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Abstract 

Motivated by the observed coordination of nearby beating cilia, we use a scale model experiment 
to show that hydrodynamic interactions can cause synchronization between rotating paddles driven 
at constant torque in a very viscous fluid. Synchronization is only observed when the shafts 
supporting the paddles have some flexibility. The phase difference in the synchronized state depends 
on the symmetry of the paddles. We use the method of regularized stokeslets to model the paddles 
and find excellent agreement with the experimental observations. We also use a simple analytic 
theory based on far-field approximations to derive scaling laws for the synchronization time as a 
function of paddle separation. 
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I. INTRODUCTION 



One of the central aims in the field of cell motility is to understand how a collection of 
beating cilia coordinates, or, on a larger scale, how a collection of swimming organisms form 
coherent patterns. For example, Paramecium swims by propagating waves of ciliary beating 
along its surface [1]. The alga Chlamydomonas beats its two flagella in synchrony to swim 
straight and asynchronously to change its orientation . At the level of a population of 
cells, sea urchin spermatozoa spontaneously form vortex patterns in the absence of cell sig- 
naling 4J . Coordination of cilia is also important in the transport of fluid. The coordination 
of nodal cilia in developing vertebrate embryos has been implicated in the determination of 
left-right asymmetry of the organism [5]. The cilia lining the human airway must beat in a 
coordinated manner to sweep foreign particles up the airway. Beating cilia may also play a 
role in the transport of sperm and egg during fertilization in mammals |(|. 

These examples are instances of the general tendency for the emergence of synchronization 
in a broad array of physical and biological systems In this article we investigate the long- 
standing hypothesis that the coordination observed in nearby beating cilia or swimmers is 



due to hydrodynamic interactions between these objects |9_|. In recent years 



;here 



been many computational and theoretical studies to support this hypothesis 
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191 ] . The key physical fact underlying all of these studies is that at the 



small scale of the cell, where the Reynolds number Re < 1, the velocity field arising from 
a deforming body falls off slowly with distance, leading to significant hydrodynamic forces 
between nearby bodies. Furthermore, the development of a fixed phase difference between 
two bodies — phase-locking — requires some kind of compliance in which the deforming body 
can adjust its beat pattern in response to hydrodynamic forces from other nearby bodies. 

The nature of this compliance is subtle. In the case of two rotating rigid helices driven 
with fixed torques (the "deformation" here is rotation), the freedom of the phase of each helix 
to speed up or slow down to maintain the fixed torque for all phase differences does not lead 
to phase-locking j^. Theoretical calculations suggest that additional degrees of freedom 
are required for phase-locking, or synchronization. For example, synchronization develops if 
the shafts of the rotating helices are connected to fixed points by stiff springs, allowing the 
axes of the helices to translate or tilt [141 ] . The directions of these small motions depend 
on whether the hydrodynamic forces are attractive or repulsive, which in turn depends in 
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FIG. 1: Schematic of the model system for hydrodynamic synchronization. Left: A pair of sym- 
metric paddles in a fluid with viscosity ry are rotated with constant torques Mi and M2 ■ The shafts 
are rigid but have flexible couplings that allow the paddles to tilt. Right: asymmetric paddles. 



detail on the phase difference (cf. the case of nearby swimmers 2l|, l22j]). 

The complexities of designing experiments that include both hydrodynamic interactions 
and controlled elastic deformation at very low Reynolds numbers have hindered experimen- 
tal studies of hydrodynamic synchronization; therefore, we built a scale model system that 
captures the essential physics, allows for detailed measurements, and is amenable to model- 
ing. This article presents results from experiments ( ffTT|) . numerical simulations ( §IIIj) . and 
a theoretical model ( §IV l and fV]) that together outline a coherent framework for describing 
hydrodynamic synchronization. 



II. EXPERIMENT 



Figure [T] illustrates the experimental configuration. Two thin paddles are immersed 
in a large tank (60 x 60 x 60cm) filled with a viscous fluid (77 = 110Ns/m 2 ), separated 
at their closest approach by a small gap, 5 = 3.6 mm. We study two different paddle 
configurations: symmetric and asymmetric. The symmetric paddle has the axis of rotation 
through the paddle center and dimensions h = 60 mm, w = 2R = 30 mm, and thickness 
t = 6 mm. The asymmetric paddle has the axis of rotation through one edge and dimensions 
of 60 x 20 x 6 mm. The paddles are small compared to the size of the tank. By repeating some 
of the experiments with the paddles at different positions within the tank, we confirmed that 
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the side walls did not affect the results in any appreciable manner. 

The paddles are supported by shafts that are hardened steel, of diameter 6.35 mm and 
length I = 120 mm, connected to the motors via flexible couplings that allow the paddles 
to tilt. The shafts are so rigid that bending due to hydrodynamic forces is negligible, but 
the couplings act as torsional springs with spring constant = 8000 mN-m/rad, leading 
to an equivalent spring constant for lateral shifts of the paddles of k = k^/i 2 - This flexi- 
bility allows the paddles to tilt slightly in response to hydrodynamic forces. We also tested 
shafts without an intermediate coupling, in which the ability of the paddles to tilt effec- 
tively vanished. The bearing assemblies are supported on separate stages to minimize any 
mechanical communication beyond hydrodynamic interactions j3] , and to allow for precise 
control of the distance of closest approach, 5. Since 5/h <C 1, the resultant flow is mostly 
two-dimensional, in the plane perpendicular to the axes of rotation. 

The two paddles are driven at constant torque using a DC servo motor, digital encoder, 
load cell, and feedback controller. Each paddle is driven by a servo motor which is encased 
in a housing. To measure the torque delivered by the motor, the housing is supported by 
bearings and prevented from rotating by a rigid, ~ 10 cm-long torque arm. Due to the 
bearings, the entire reaction torque on the housing is transmitted by the torque arm to a 
precision load cell. The load cell output signal is used as a feedback to a PID controller that 
adjusts the voltage driving the servo motor, thus maintaining a defined torque. The PID 
controller updates at approximately 100 Hz — 500 times faster than the typical rotational 
frequency of the paddles in the experiment (0.2 Hz). The position of the paddle is recorded 
from the output of the digital encoder at each update of the PID controller. Velocity is 
calculated from the position using high-order finite differences. The system was calibrated 
by measuring the rotational speed vs. voltage for an isolated paddle over a range of torques 
and using the theoretically-known torque-speed relationship to associate the measured load 
cell voltage with a specific torque. The accuracy and stability of the system was verified 
by measuring (i) the torque fluctuation for an isolated paddle rotating at constant speed, 
and (ii) the velocity fluctuation of an isolated paddle rotating at constant torque. In both 
configurations, we confirmed that the system was stable to better than 1.5% of the set 
point. Typical driving torques range from 4mN-m to 25mN-m, corresponding to rotation 
frequencies no more than 0.2 Hz. 

At these conditions, the Reynolds number, Re = pu>R 2 /r] ~ 10 -3 , is small enough to 
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FIG. 2: (Color online.) Phase difference A9 = 6 2 — 0\ vs. dimensionless time t/To (Tq = Qttt]E? / M) 
for symmetric paddles with (a) Mi = M 2 and stiff shafts, (b) (M 2 - Mi) /M 1 r* 0.003 and stiff 
shafts, and (c) M\ = M 2 and shafts with flexible couplers. 
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FIG. 3: (Color online.) Dimensionless synchronization time vs. dimensionless gap size 5/R for 
symmetric paddles (circles) and asymmetric paddles (squares). For symmetric paddles, the dimen- 
sionless synchronization time is measured from the moment of phase difference A9 = 0.1 to the 
time of the first stable state A9 = ir/2. For asymmetric paddles, the dimensionless synchroniza- 
tion time is defined as the time from A8 = 0.8 to A9 = 0. The uncertainty comes from the low 
frequency fluctuation in AO due to the system noise. Time is normalized by To = QirriE? / M . 

justify the neglect of inertial forces. This was confirmed experimentally by noting that if 
the paddle rotation was initiated with a constant velocity, the time taken to reach constant 
torque was less than 250 ms. For these Stokes flows, the characteristic velocities scale linearly 
with the motor torques (Mi,M 2 ), and the state of the system is determined by the angles 
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FIG. 4: (Color online.) Phase difference vs. dimensionless time for symmetric (red line) and 
asymmetric paddles (green line), compared to simulation (blue lines). Time is measured in units 
of To = 67rr]R 3 /M . The insets show the phase difference once phase-locking is achieved. In both 
cases, the normalized gap between the paddles is 5/R = 0.24. 

0.5 

0.4 

-5.0.3 
as 

< 

0.2 
0.1 

°'°0 200 400 600 800 

t/T 

FIG. 5: (Color online.) Phase difference vs. dimensionless time t/To for symmetric paddles driven 
at torques M 1 = M 2 = 8.4, 12.8, 17.5, and 22mN-m, with 5/R = 0.24. The collapse of the data 
shows that the time to synchronize scales with To = 6irr]R 3 /M. The inset shows the dimensionless 
period T/Tq in the synchronized state. 

of the two paddles {9i,9 2 ) (Fig. [Q and the small shifts of the paddles due to the flexible 
couplers. In the high-stiffness case, the paddles did not synchronize in any measurable 
time; instead, the phase of each paddle increased roughly linearly with driving torque, 
(0i,#2) oc (tMi,£M 2 ), independent of the initial phase difference (Fig. [2|l, b). However, 
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paddles with flexible couplers and small distance of closest approach locked phases in 10-20 
revolutions (Fig. [2b). Note that we measure time in units of To = 67T1]R 3 /M, (M is the mean 
torque) which is roughly one tenth of a rotation period. The data we display in this article 
is for a dimensionless gap size 8/R = 0.24. We also varied 5/R for both kinds of paddles 
from pa 0.1 to ~ 0.6, and found that the time to synchronize increased with spacing, with 
longer times and a faster increase for the symmetric paddles (Fig. [3]). 

For Mi = M2, the symmetric paddles locked phases at A9 = 82 — 61 = 7r/2, and the 
asymmetric paddles settled at A9 = (Fig. H]). These two states represent the conditions 
that roughly maximize the distance of closest approach of the two paddles. Since the paddles 
would minimize their distance of closest approach if they maintained their typical initial 
phase differences (A9 = for the symmetric paddles, A9 = ir for the asymmetric paddles), 
the rotation speed of each paddle rises as the paddles synchronize. Denoting the rotation 
speed of an isolated paddle by u , we found that the speed of both symmetric paddles rises 
from 0.72uq to 0.85c<Jo as synchronization develops, whereas the speed of both asymmetric 
paddles rises from 0.75co>o to 0.93co> . While these synchronized states are stable, there is 
a consistent and repeatable phase fluctuation (Fig. Hj-inset) corresponding to the variation 
in rotational speeds as the hydrodynamic interactions between the paddles wax and wane 
during a cycle. The fluctuation amplitude in the asymmetric case is larger than in the 
symmetric case because there is a larger variation in the distance between the asymmetric 
paddles during a period. These observations qualitatively agree with the results of numerical 
calculations on rotating rigid helices with flexible couplers fl4j |. In our experiments the 
phase fluctuations and rise in velocity as synchronization develops are more dramatic since 
the variation in the hydrodynamic interaction between paddles over a period is greater than 
in the case of helices. 

The final state of synchronization was found to be independent of the initial orientation 
of the paddles. The time to synchronize scales with To, perhaps with a weak dependence 
on torque (Fig. [5]). The number of paddle revolutions needed to synchronize is therefore 
roughly constant, 15 in the case of symmetric paddles, and 20 for asymmetric paddles. In the 
synchronized state, however, the dimensionless rotation period T/T increases slightly with 
torque (Fig. [5j inset). When the symmetric paddles are operated with a torque mismatch 
between the two motors, the synchronized phase difference increases with AM, although for 
a large mismatch, AM /Mi = 3%, the synchronized state is only marginally stable and the 
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FIG. 6: (Color online.) Phase difference vs. dimensionless time for symmetric paddles with 
5/R = 0.24 and AM /Mi = (M 2 - M 1 )/M 1 = -1% (a), AM /Mi = 0% (b), AM /Mi = 1% (c), 
and AM /Mi = 2% (d). When AM /Mi ~ 3%, the synchronized states become unstable, and the 
system exhibits transition in which AO jumps by ir. 

phase difference can jump abruptly by A8 = n (Fig. EJ). 



III. NUMERICAL SIMULATIONS 



These experiments give strong evidence that the phase-locking of the paddles is due to 
hydrodynamic interactions. We tested this hypothesis by using the method of regularized 



Stokeslets [24J to model the flows induced by the paddles. Each paddle is replaced by a 
rectangular array of regularized Stokeslets S b with strength f u , where /i and v label the 
Cartesian coordinates x, y, and z. The flow from the Stokeslet at x' is given by 

*V(x) = ]T^(x,x')/„ (1) 



with associated pressure 



The Stokeslet S b u satisfies 



P(x) = ^^^- (2) 



Y,dSlJdx v = (3) 

V 

V 2 ^(x, x') - dpi/fa^ = -8tt V0 6 (x - x'), (4) 
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where 06 (x — x') is a smooth approximation to the Dirac delta function with spread b, 



and r = |x— x'|. The number of stokeslets and the spread b are chosen to give good agreement 
between the measured and simulated resistance coefficient for a single rotating paddle at 
the center of the tank. The spread b is large enough to make the regularized stokeslets 
overlap, which prevents fluid from leaking through the paddles. We model the flexibility of 
the couplers with springs of torsional spring constant fcj- For simplicity we suppose that 
the shafts are always vertical, but can undergo slight shifts in the horizontal plane. With 
the assumption that the paddles are rigid, the degrees of freedom are the angles (#1,6*2) of 
the paddles and the positions of the shafts. Balancing forces and torques leads to coupled 
nonlinear differential equations which we solve numerically. Figure H] shows the excellent 
agreement between the experiments and the simulations for both the asymmetric and the 
symmetric paddles. The simulation accurately captures the frequency and amplitude of the 
oscillations associated with the rotation of the motors, as well as the slower evolution of the 
phase-locking. When the driving torque is varied over the range used in the experiment, the 
simulations yield that the dimensionless time to synchronize T s /Tq remains approximately 
constant, with a weak dependence on torque, in accord with Fig. [51 Simulations with infinite 
spring constant &t show no phase-locking. Since the paddles in the simulation are coupled 
only through the hydrodynamic interaction, we conclude that the cause of the phase-locking 
is the hydrodynamic interaction and not any stray mechanical coupling that might be present 
in the experimental apparatus. 

IV. SIMPLE MODEL FOR ASYMMETRIC PADDLES 
A. Oseen tensor model 

We can gain more insight into the mechanism of phase-locking by developing a simple 
theory along the lines of reference [3]. A minimal model for the asymmetric paddles is 
to replace each paddle with a sphere of radius a attached to one end of a rod of length R 
(Fig. [7]). The rod is rigid and does not disturb the fluid. The other end of the rod is attached 
to a stationary point by a spring with spring constant k. The rods are rotated by moments 
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FIG. 7: Model for asymmetric paddles. The figure is not to scale; note that o<i?<0. 



Mi and M 2 which are applied at the ends of the rods attached to the springs, where we can 
imagine shafts perpendicular to the plane of the page. The spring is stiff, with k 3> Mi/ R 2 . 
Denote the positions of the balls by = =f(.D/2)x + x, + Rpi, where the minus sign applies 
for i = 1, the plus sign applies for i = 2, and pi = (cos0j,sin0j). Note that Q{ is defined as 
the angle the rod makes with the x-axis, not the angle rj makes with the x-axis. The vectors 
xx = (xi,yi) and x 2 = (x 2 , 2/2) are the displacements of the shafts from the stationary points 
(—D/2,0) and (D/2,0), respectively. If we suppose the balls are far apart, with D ^> a, 
then the leading-order interaction between the two balls is given by the Oseen tensor 251 ]: 
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677770, 
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(f 2 • r 12)1*12 

I r 12 1 3 
(fi ■ ri 2 )r 12 

I r 12 1 3 



(6) 
(7) 



where v, = dr^/dt = is the velocity of the ith ball, fi and f 2 are the forces exerted by the 
balls on the fluid, and ri 2 = ri — r 2 . 

Since the spring is assumed linear and the motion of the rod incurs no drag force, the 
balance of forces on each paddle is — f, — A;Xj = 0. We must also enforce moment balance. 
Since inertia is unimportant at Re = 0, we may compute moments about the points Xj for 
each paddle: 

Mi + z • {Rpi) x (-ft) = 0. (8) 



B. Separation of time scales 

The analysis of the equations of motion is simplified by the recognition that our problem 
has three well-separated time scales: (i) a short time scale = rja/k (recall k = k^/i 2 ) 
that controls the rate of relaxation of the springs, (ii) an intermediate time scale T\ = 
QurjaR 2 / Mi [28| that controls the period of an isolated ball driven by torque Mi, and (iii) a 
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long time scale T s that characterizes the time for phase-locking to develop. Since the springs 
are stiff, Tk <C T\. Since the interaction between the paddles is weak, T\ <C T s . The goal 
of the simple model is to give a clear derivation of how T s depends on the parameters of 
the problem. Since the phase-locking arises from hydrodynamic interactions, which vanish 
when a/D — > 0, we expect T s to scale as some power of D/a for large D/a. 



C. Dimensionless far- field equations of motion for 9i and Xj 



Since the balls are far apart, we expand the equations of motion in powers of 1/D, 
assuming that a ^ D and R <C D. Measuring length in units of R, time in units of T\, and 
using force balance to eliminate fj leads to the dimensionless equations of motion, 

(X + 0) = -HX/e. (9) 

In Eq. USD, e = M l /(kR 2 ) < 1, X and are 4x1 vectors with 



X 





|. *-| 
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(10) 



where 0, 



sin^j, cos^j). The 4x4 matrix H is the Oseen tensor to leading order in a/D, 
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l + xx) 
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(11) 



where I is the 2x2 identity matrix and xx is the 2x2 matrix with unity in the upper left-hand 
corner and zeros elsewhere. 

Using Mi as the unit for torque, the moment balance equations (JSJ) take the form 

1 + (xi/e) • 0i = (12) 

(13) 



AM 

1 + — + (x 2 /e) • 9 2 = 0, 



Mi 

where AM = M 2 — Mi. From these equations we conclude that Xj is 0(e). Note that the 
shafts have a nonzero displacement Xj even when the paddles are isolated. 



D. Far-field equations of motion for average angular speed and phase difference 

To understand phase-locking, it is not necessary to resolve the motion of the paddles on 
the short time scale T^. In dimensionless variables, these short-scale motions are character- 
ized by transients of the form exp(— t/e). By considering dimensionless times i>ewe may 
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neglect these transients and treat eX as small. Physically, this approximation reflects the 
fact that once the transients have decayed, the drag forces incurred by the small motions X 
arising from the extension of the springs are small, but not negligible, compared to the drag 
forces due to the rotation © of the balls about the shafts. Therefore, we solve Eq. (Q for X 
using iteration, finding 

X^ -eH^e + e^- 1 ^ (V 1 ©) . (14) 

In terms of i: we have 

X » e (~ 0A + ~^ A \ + e 2 [~ bVpX + (15) 

where we have only retained terms of 0(a/D). 

In Eq. (1151) . we have discarded terms of the form Oi, since they are 0(a 2 /D 2 ). To see 
why, observe that for time scales longer than the motion is characterized by two well- 
separated time scales, T\ and T s . The form of the interaction suggests that T s oc D/a. To 
explicitly account for the multiple scales Ti and T s , write [26[] 

1)2 = w(t)*tA0(t)/2, (16) 

where r = at/D describes the slowly- varying time dependence of the rotational frequency 
and the phase difference. Note that uj(t) is the average angular speed, and AO is the average 
phase difference. The angular speed and phase difference also have rapidly vary parts with 



zero average, but these are lower order in a/D 261 ] . Equation ffl6l) shows that the leading 
term of 9i is (ci/D)uj'(t). But since the average rotation speed to is constant in the absence 
of interactions, uj'(t) must be at least 0(a/D). Thus, 6*j is at least 0(a 2 /D 2 ). 

To find the governing equations for angular speed to and phase difference AO, substitute 
the shaft displacements x» from Eq. ( fl5l) into moment balance, Eqs. ( Tl2lTl~3]) . Finally, average 
the resulting equations over a period, treating the slowly-varying variables u and AO as 
constants under the average. We find that the average dimensionless speed is given by 

AM 9 Ou 

w = mT + I7) cosM ' (17) 

The interacting paddles turn faster than they would in isolation. This result is in contrast 
with our paddle experiments, where we saw in ^JTT] that the asymmetric paddles rotated more 
slowly compared to an isolated paddle. It is too much to demand that our far-field theory 
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captures every aspect of the paddle experiments, since the paddles are close to each other 
in the experiment and the theory is valid when they are far apart. 
The dimensionless phase difference obeys 

dA# 9 a . A „ AM 

— — = — e— sinAfl + — — . 18 
dt 2D Ml k j 



These results flTTHTHJ) are equivalent to the results of reference [18j. For equal driving torques, 
AM = 0, Eq. ( TlBl shows that the paddles synchronize to A9 = 0, independent of the initial 
value of A9, in (dimensional) time T s ~ (D /a){kR 2 /M\)T\, or 

^ DkR HiryaR 2 
s ~ a Mi Mi 1 ' 

When Mi ^ M2, the paddles phase-lock with a nonzero phase difference, which increases 
to 7r/2 in the steady state as the torque difference increases to the critical value given 
by AM /Mi = (9/i)(a/D)M 1 /(kR 2 ). Note that the factor of a/D and the smallness of 
e = Mi/(kR 2 ) mean that M 2 must be very close to Mi for the phase difference A9 to have 
a fixed point. Thus, in the derivation of Eqs. ( JT7HT8I) we considered AM /Mi and TiAO to 
be 0(ea/D). 

This simple theory predicts that T B /T\ varies inversely with torque, whereas the exper- 
iments show that T s /Ti depends at most weakly on torque (Fig. [5]). Again, the resolution 
of this discrepancy is that the simple theory is valid in the far-field limit with D ^> a, 
whereas the experiments are carried out in the near-field regime where T s /Ti is independent 
of torque. 



E. Physical explanation for phase locking 

Each of the terms of Eq. ( fl5l) has a simple interpretation. First consider the limit of an 
isolated paddle, a/D = 0. To leading order in e, the ball on the end of the rod undergoes 
circular motion. This motion leads to a drag in the —9i direction, which stretches the spring 
along —9i, which in turn leads to an 0(e) component of the ball's velocity parallel to the 
rod, along the direction (see the left ball in Fig. [Hh). In our dimensionless units, the ball 
exerts an 0(1) force on the liquid in the 9i direction, and an 0(e) force on the liquid in 
the pi direction. To get the displacement X, we multiply these forces by e, and thus get 
the 0[(a/D) ] terms of Eq. ( fl5i) . Now consider the hydrodynamic interactions. For a given 
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FIG. 8: (Color online.) Physical explanation for synchronization. The figure is not to scale and 
the paddles have been moved artificially close, (a) On the left we show the components of velocity 
of the first ball (blue arrows). On the right, we show the components of the drag induced on the 
second ball by the motion of the first ball (red dotted arrows), and the components that contribute 
to the hydrodynamic torque on the second paddle (green arrows), (b) The same situation as (a), 
but showing the velocity components of the second ball and the induced forces on the first. The 
difference in the hydrodynamic torques tends to make A8 = 0. 

paddle, each of the forces just described induces a Stokeslet flow, falling off inversely with 
distance, and leading to the 0(a/D) terms in Eq. (IT51) . 

The drag forces on each paddle induced by the motion of the other are shown in Fig. [HJ 
From this figure we can see why the paddles synchronize. Suppose that the second paddle 
slightly lags the first. Since the spring is flexible, the ball of the paddle on the left has a 
velocity component of ed\ along the rod as well as the component Q\ perpendicular to the 
rod (blue arrows, Fig. [Ht, left). This motion induces drag forces on the ball on the right (red 
dotted arrows, Fig. [HJi, right), which in our dimensionless units are down by a factor of a/D 
from the velocities. The components of these forces perpendicular to the rod (green arrows, 
Fig. [HK, right) contribute to the hydrodynamic torque on the paddle. Likewise the motion of 
the paddle on the right (blue arrows, Fig. Eb, left) induces forces that lead to hydrodynamic 
torques on the left paddle. The phase difference A8 is governed by the difference of the 
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FIG. 9: Model for symmetric paddles. The figure is not to scale; note that a <C R <C D. 

torques, which for small A6 is given by the difference of the small (green) arrows in Figs. 
and b. The torque difference makes A6 = a stable fixed point (for AM = 0). 

F. Power dissipation 

We may readily examine the question of power dissipation using our simple model. First 
note that for fixed driving torques, the power dissipated decreases when the hydrodynamic 
resistance of the paddles increases. Therefore, when AO = 0, the drag is minimized and 
the dissipation rate is maximized. As AM increases, the increase in AO leads to greater 
resistance and therefore lower dissipation rate. To leading order in e, we may use Eq. ( TPTj) 
to show that the dimensionless power averaged over one period, P = M\0\ + M262, takes 



In general, the phase difference chosen by the system does not minimize the power dissi- 
pated. The same conclusion has been reached for the hydrodynamic phase-locking of nearby 
swimming sheets 

V. SIMPLE MODEL FOR SYMMETRIC PADDLES 
A. Oseen model and nondimensionalization 

To understand why the symmetric paddles lock phases with AO = tt/2 when AM = 0, 
we model the paddles as dumbells (Fig. [9]). Each dumbell consists of two balls connected 
by a rod that does not disturb the fluid as it moves. The midpoint of each rod is attached 



the form 




(20) 
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to a fixed point by a stiff spring, and the balls at the ends of the rods have positions 



r 1± = -{D/2)± + x. 1 ±R Pl (21) 
r 2± = (D/2)x + x 2 ±Rp 2 , (22) 

where x 4 is the displacement of the midpoint of the ith rod from the corresponding fixed 
point. Denoting by fj± the forces that the balls on the zth dumbell exert on the fluid, the 
balance of forces on each dumbell implies 

- f i+ - f-_ - k Xi = 0, (23) 

and the balance of torques implies 

Mi + z • (Rpi) x (-f i+ ) + z ■ {-Rpi) x (-$_) = 0. (24) 

Assuming all balls are far apart, we again use the Oseen model, Eq. ([7j), this time extended 
to the four balls labeled a = 1—, 1+, 2—, 2+: 



-E 



67rr)a 8711] L 



fa/3 1 | r a/3 
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(25) 



where r a/ g = r a — rp. This is valid when R 3> a and D ^> a, but we will also assume D R. 

The 0j I— > 0j + 7r symmetry of the dumbells makes the hydrodynamic interaction between 
the dumbells more subtle than the asymmetric case. First observe that the spring of an 
isolated rotating dumbell does not stretch since the net hydrodynamic force on the balls 
vanishes. Thus, Xj = when D — > 00. However, for finite D/R, the flow induced by the 
rotation of one dumbell causes the spring of the other dumbell to stretch. To estimate the 
amount of stretch, consider the flow induced by dumbell 1 at dumbell 2. The far-field flow 
is an asymmetric force dipole, also known as a rotlet, falling off inversely with the square of 
distance |2?J] . Thus, the flow v 2 % induced at dumbell 2 is approximately v 2 i ~ f\R/ (r]D 2 ), 
leading to drag on dumbell 2 of about rjav2i ~ af\R/D 2 ~ aM\/D 2 . This drag causes the 
spring of dumbell to stretch, with a displacement 

|x 2 [ ^R^M^ 

R ~ RD 2 kR 2 ' 1 ' 

As in the previous section, it is convenient to measure length in units of R and time in 
units of T\ = 6nr]R 2 a/Mi. Thus, the dimensionless displacement is |x 2 | ~ ea/D 2 . 
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B. Far-field equations of motion 



For each spring of a pair of rotating dumbells, the leading order stretch of the spring is 
second order in D^ 1 . However, to derive equations describing phase-locking of symmetric 
dumbells, we will see that we must expand the displacements x, to 0(D~ 3 ). These third 
order displacements arise from reflections of the dipole force. For example, the (D(D~ 2 ) 
deflection of spring 2 from the dipole originating at dumbell 1 induces a point force at 
dumbell 2. This point force causes an 0(D~ 3 ) displacement at dumbell 1. As we will see, 
to find x,j to 0(D~ 3 ) we need only expand the Oseen tensor to second order in D~ x . Thus, 
to second order in -D -1 , Eq. (|2"5|) becomes 



Xl - 0101 
±2 + O 2 02 

\x 2 - e 2 o 2 j 



H 



/f 1+ \ 
fi 
f 2+ 
VW 



(27) 



where H is the 8x8 Oseen tensor evaluated at Xi = x 2 = and expanded to 0(D~ 2 ). 

Before describing H, it is convenient to re-express Eq. ( 1271) in terms of the sums and 
differences of forces on each paddle, f, = fj + + fj_ and Afj = f i+ — fj_, respectively. Note 
that in terms of these variables, the dimensionless force-balance equation (|23|) becomes 

- fi - Xi/e = 0, (28) 

and the dimensionless moment-balance equations become 

(29) 
I ^— - Af 2 • § 2 = 0, (30) 



i - Afi ■ e l 

AM 



Mi 

where Mi is the unit for torque. Returning to Eq. ([2~7]) . we add and subtract the appropriate 
rows of Eq. ( 1271) and rearrange to find 



X = -AX + BAF 

e 

= -CX+DAF. 

e 



(31) 
(32) 



where X and are defined as before in Eq. (TlOl) ; the 4x4 matrices A, B, C, and D are 
given in the appendix; and 

w 

Af 2 



AF 



(33) 
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Expanding in powers of D l , we find (see Appendix) 



A = A (0) + D _1 A (1) + 0(/r 3 ) (34) 

B = D~ 2 B^ + 0(D- 3 ) (35) 

C = ZT 2 C (2) + 0(D- 3 ) (36) 

D = D (0) + 0(/r 3 ). (37) 

Likewise, we expand X and AF in powers of D^ 1 : 

X = X (0) + D^X® + ZT 2 X (2) + • • • (38) 
AF = AF (0) + D^AFW + £)- 2 AF (2) + ■ • • . (39) 

Since B and C are 0(D~ 2 ) at leading order, our order of expansion is sufficient for deter- 
mining X to (D(D~ 3 ) and AF to 0(D~ 5 ). At zeroth order, we find X^ = 0, as expected, 
and AF<°> = [D^]- 1 ^, with 

Aff } = (40) 

l-3a/8 V ; 

Substituting Af(°) into the moment balance equations (129"H3"U1) . taking their sum and differ- 
ence, using Eq. f fT6|) to eliminate 9\ and 82 in favor of u and AO, and integrating over one 
period yields the dimensionless average speed and phase difference, 



UJ 



(0) 



l/2(l-3a/8) (41) 
A# (0) = (l-3a/8)AM/Mi. (42) 

Since lu^ and AQ^ are independent of D, they are the average speed and phase- difference, 
respectively, for non-interacting dumbells. There is no phase-locking if there is no interac- 
tion, and the phase difference increases in proportion to the difference in driving torques, 
AM. 

Note that the factors of 3a/8 in Eq. (140]) are due to the interaction between the two balls of 
a given dumbell: one ball induces a disturbance flow of magnitude (67rrja9) / (87rr/2) = 3a#/8 
at the other ball. This disturbance flow hinders the motion of the other ball. 

The leading order displacements of the shafts are given by 

X (2) = i A (0) X (2) + B (2) AF (0)_ (43) 
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As in the case of the asymmetric paddles, this equation is readily solved to 0(e 2 ); however, 
the full expression is so cumbersome that we only report the result to leading order in e and 
a in the appendix. The next order contribution to the force difference is given by 

= I C (2) X (2) + D (o) AF w (44) 
e 

Again, the full expression for AF^ 4 ) is so cumbersome that we only report the leading order 
terms in the appendix. Using AF^ 4 ) in the difference of the moment equations and averaging 
yields terms proportional to AO, which do not lead to phase- locking. The average of the sum 
of the moment equations leads to a decrease in the average rotation speed, which together 
with Eqn. (jH]l yields 

1 3a 153 a 2 
" = 2-16-l6^- (45) 

The interacting paddle turn more slowly than they would in isolation. 
The third-order displacement of the shafts is determined by 

X(3) = Ia^X^ + -A(°'X( 3 '. (46) 
e e 

Solving for X^ 3 ) (see Appendix for leading terms), and substituting into 

= ~C^X® + D^AF® (47) 

yields AF^ 5 ) (see Appendix for leading terms), which has terms that lead to phase-locking. 
Using moment balance Eqs. (129H301) and averaging, together with the leading order result 
(SOD, yields 

1 243 a 3 

AO = - AM + — e— sin 2 AO (48) 

2 8 D 

Equation fj4*8l is the main result of this section. The (dimensional) synchronization time 
for the symmetric paddles scales as 

-D 5 kR 2 67T7]aR 2 

When AM = 0, Eq. (T481) has a stable fixed point at AO = ir/2, in accord with our experi- 
ments and the more accurate regularized stokeslet simulation of §1111 As in the case of the 
asymmetric paddles, the torque difference must be small for phase-locking to occur. The crit- 
ical torque difference, above which phase-locking cannot occur, is AM crit = (243/4)ea 3 /D 5 . 
Note that the average phase-difference in the phase-locked state depends on AM. Note also 
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that the time for phase-locking depends more strongly on separation for the symmetric pad- 
dles compared to the asymmetric paddles. It is not easy to give a simple physical picture for 
why the paddle separation D enters the synchronization time with a fifth power. We simply 
note two effects: (1) the flow induced by the force dipole of one paddle reflects off the other 
paddle, and then again off the first paddle, leading to four powers of -D -1 , and (2) the torque 
exerted by a flow on the paddle arises from the difference in the flow at the two ends of the 
paddle, leading to another factor of D^ 1 . Our case is reminiscent of the fifth power that 
appears in the reorientation of oscillating dumbells 2l|. Although our experiments were not 
carried out in the far-field regime, we found that the synchronization time depends more 
strongly on separation in the symmetric case compared to the asymmetric case (Fig. [3D- 
Finally, we note that to leading order in D^ 1 , the power dissipated in the synchronized 
state is independent of A8, since the average rotation speed u in the synchronized state is 
independent of A8 (Eqn. 145]). 



VI. CONCLUSION 



To summarize, we have presented perhaps the simplest experimental realization of the 
phenomenon of hydrodynamic synchronization at low Reynolds number. The requirements 
for synchronization are subtle: the system must have a slight flexibility to allow small shifts 
in the positions of the paddles. Since this flexibility is generic, we expect that conditions 
allowing hydrodynamic synchronization will commonly arise in a wide range of systems at 
low Reynolds number. On the other hand, our work indicates that hydrodynamic synchro- 
nization is not robust, since it requires that the driving moments be fine-tuned to be close 
to each other. 
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APPENDIX A: OSEEN TENSOR FOR SYMMETRIC PADDLES 



For H evaluated at x» = and the matrices A, B, C, and D defined in Eqs. (I3TH32I) . 
expanding in powers of D~ x yields 



(-\ 



32 

— ■§2 a sin(2#i 



a cos(2#! 



3a 
'2D 



-^asin(2#i) 
-§§ + ^acos(20i; 


_ 3a 
AD 



3a 
2D 



\ 





ff - ^acos(20 2 



3a 
'AD 



32 32 

-4asin(20 2 ) 



~asin(20 2 ) 



|f + ^acos(20 2 )/ 



and 



D 



9o 
32 



— -^a cos(20i 
^asm(28 1 ) 



32 



a sin(26?i 



V 








i_| + | acos ( 2 1 










If - |acos(20 2 ) 

_3 

32 



-:|asin(26» 2 ) 



(Al) 




-^asin(20 2 

9a , 3 



I_| + A acos( 20 2 )/ 
(A2) 

for the blocks on the diagonal of the Oseen tensor. For the blocks off the diagonal, we have 



B 













-11? cose 

3a 2 
AD 2 



3a 



in 0o \ 



3a 2 
2D 2 



sin 2 COS 2 



C 



f^cos0! -|^sin0! 
V-gsin0 1 |^cos0 1 

/ 


Wi sin 2 

3a cos 0o 








—Ms cos 0i 








(A3) 



3a 
AD 2 



sin0! ^ 



2^2 COS 2 

3a sin 02 



2D 2 

^sin0! cos 0i 







(A4) 



AD 2 ° 2 AD 

The second order spring deflection to leading order in e and a is 

/ (9/2)02 sin(20 2 ) \ 
-30 2 



X (2) 



(If 



-(9/2)0! sin(20!) 
30i 



0{ta z 



(A5) 
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The fourth order force difference to leading order in a and e is 

/-0i [5 sintfx + 3sin(30i)]\ 

A , 41 9 2 29 1 cos9 1 (1 + 3 sin 2 0j 
AF (4) « -a 2 v y 

4 -02[5sin0 2 + 3sm(30 2 )] 
\ 20 2 cos 2 (1 + 3 sin 2 2 ) / 

The third order spring deflection, leading order in e and a, 

[ 30isin(20i) \ 

9 



X (3) 



-a 2 e 



\ 



-0i 
-30 2 sin(20 2 ) 

2 



Finally, again to leading order in e and a, 

[-02 [sin 0x + 6 cos 0i sin(20 2 )]\ 

27 , 



AF (5) 



2 [cos 0i + 3sin0! sin(20 2 )] 
-0i [sin 2 + 6 cos 2 sin(20i) 
\ 0i [cos 2 + 3 sin 2 sin(20i)] ) 



(A6) 



(A7) 



(A8) 
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